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ABSTRACT 

We estimate a value of the viscosity parameter in AGN accretion discs for the 
PG quasar sample. We assume that optical variability on time-scales of months 
to years is caused by local instabilities in the inner accretion disc. Comparing the 
observed variability time-scales to the thermal time-scales of a-disc models we 
obtain constraints on the viscosity parameter for the sample. We find that, at a 
given L/LEdd, the entire sample is consistent with a single value of the viscosity 
parameter, a. We obtain constraints of 0.01 < a < 0.03 for 0.01 < L/L^dd < 1.0. 
This narrow range suggests that these AGN are all seen in a single state, with 
a correspondingly narrow spread of black hole masses or accretion rates. The 
value of a we derive is consistent with predictions by current simulations in 
which MHD turbulence is the primary viscosity mechanism. 
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1 INTRODUCTION 

The transport of angular momentum remains a major 
unknown in accretion disc theory. Shakura & Sunyaev 
(1973) suggested that magnetic fields and turbulent mo- 
tions are likely candidates for the angular momentum 
transport. They argued that turbulence is a dominant 
factor and postulated that the stress is proportional 
to the local pressure. The proportionality parameter a, 
which parametrizes the underlying transport process, de- 
scribes the efficiency of the turbulent transport. The a- 
parameter also represents the scale of the largest possible 
turbulent cell. The a parameterization provided a base 
for developing the accretion disc theory and linked the 
theory with the observations. 

More recent work has shown that in weakly mag- 
netized discs MHD turbulence can provide the major- 
ity of the outward tr ansport of angular momentum 
fealbus fc Hawlevll99J) . Balbus & Papaloizou (1999) re- 
viewed the a-disc models and found that MHD turbu- 
lence follows the a prescription very closely. They show 
that the vertically averaged disc with viscosity due to 
magnetorotational instabilities is described by a simi- 
lar set of radial equations to the a-discs. In particular 
the local energy dissipation rate is determined by the 
global disc parameters as it is in the a-discs. However, 
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global 3-D numeric al simulations of accretion d iscs are 
extremely complex jMatsumoto fc Shibatalll99 , it) and so 
most work to date has been done using local 2-D shearing 
box, cylindrical and axisymmetric limit models. All these 
methods have provid ed predictions of the cv-p arameter 
ranging from —0.005 jB randen burg et a~ 1 ll995l) to -0.6 
filawlev. Gammie fc BalbusHl995l) 

In order to determine the magnitude of the viscosity 
better estimates from larger samples of all categories of 
accreting objects are needed with independent methods 
for estimating a. Since the spectrum of a disc at any par- 
ticular instant is not strongly dependent on its viscosity 
the most effective way to estimate the viscosity is through 
the variability of the disc luminosity. A number of studies 
of this kind have been done for accreting binary systems 
since the time-scales of variability are optimal for obser- 
vation, generally ranging from days down to a fraction of 
a second. 

In the case of AGN we observe typical optical fluc- 
tuation time-scales of weeks to years although there 
have been only a few long-term monitoring programs 
to provide such data (eg. Pica et al. 1988; Giveon et 
al. 1999, hereafter G99). At optical and ultraviolet fre- 
quencies the majority of the continuum luminosity seen 
in AGN is held to be em ission from the accretion disc 
iMalkan fc Sargendll982T) and large variability has been 
observed in these wav ebands for many such objects 
iWebb fc Malkanll2000l) . The physical origin of optical 
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flux variations has not been pinned down, but if disc emis- 
sion dominates the spectrum at optical wavelengths the 
variability is expected to be related to the disc instabili- 
ties. 

In this study we use a subset of the Palomar Green 
sample of AGN (G99), spanning approximately seven 
years of data, to derive limits on the viscosity parameter 
for the the sample as a whole via the method described 
in Siemiginowska & Czerny (1989, hereafter SC89). We 
assume any fluctuations in the lightcurves are due to an 
instability arising in the innermost disc regions and can 
be associated with the thermal time-scale. We do not 
identify the exact physical process responsible for this 
instability. It can be related to the Lightman-Eardley 
instabilities in the radiation pressure dominated region 
of the disc (iLightman fc Eardlevlll973) . but our method 
does not require modelling of the instability and therefore 
could be applied to any type of instability which may be 
present in the inner disc. Using a set of multi-temperature 
blackbody disc models, based on the standard a-disc as 
prescribed by Shakura & Sunyaev (1973), we calculate 
the contribution to the optical spectrum from the outer, 
stable disc and assume that there is no contribution from 
the innermost parts of the disc where the instability has 
developed. Comparing the variability time-scales of the 
sources with the thermal time-scales of the models allows 
us to constrain the Q-parameter for the sample. 



2 DATA SAMPLE 

The data are taken from the optically selected Palomar- 
Green (PG; Schmidt & Green 1983) sample. A subset of 
42 bright (-28 < M B < -22) and nearby AGN (0.06 < 
z < 0.37) in the northern hemisphere were observed at 
the Wise Observatory, Israel (G99) in 1991-1998. Ob- 
servations were made with a CCD camera mounted on 
the lm optical telescope and span 7 years with a typical 
sampling interval of 39 days - a higher, more even sam- 
pling interval than that of previous programs (eg. Pica 
et al. 1988). The wavelength range covers the Johnson B 
band, 3600-5600 A (A central = 4400 A), which is a region of 
the spectrum believed to be dominated by the disc emis- 
sion. The PG sample is statistically complete with good 
photometric accuracy over this wavelength range (~0.02 
mag. at 4400 A). G99 assume a Hubble constant Ho = 
70 kms - Mpc~ , a deceleration parameter of qo = 0.2 
and a power-law shaped continuum, F v oc v~ a where the 
spectral index a = 0.5 for the K-correction. This is a flux- 
limited sample so we expect to see selection effects as dis- 
cussed in Schmidt & Green (1983). This sample has also 
been studied at a number of other wavelengths in clud- 
ing R band (G99), radio dKellermann et al. Ill989ft. IR 
jNeueebauer et al. 1 119871) and X-rav llLaor et al. Ill997f) 
and its optical polarization properties were studied by 
Berriman et al. (1990). Measurements of the masses of 17 
PG quasars included in th is subsample have been made 
by reverberation mapping ijKasoi et al. ll200nl^ . 

From the lightcurves we derived a two-folding time- 
scale as a characteristic time-scale for each source. The 
two-folding time-scale gives the time for an objects' lumi- 
nosity to change by a factor of two (either an increase or 



a decrease) which can be used as a measure of variability 
in both the data and the models (O'Brien, Gondhalekar 
& Wilson 1988 and references therein; SC89). 
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where /max and / m i n are the maximum and minimum 
observed flux levels respectively, r s = r max — T m i n is the 
time between maximum and minimum flux levels and 
T2 is the two-folding time-scale. The lightcurves for this 
data set do not show flux variations of a factor of two, 
in fact <jb/B <34%, in which case T2 is the time over 
which a linearly extrapolated flux variation would in- 
crease or decrease the observed flux by a factor of two. 
These timescales were then corrected for the effects of 
cosmological time dilation. 

In order to calculate / m in/(/max — /min) we used the 
standard transformation from magnitude to flux, 



m min — m max = —2.5 log 
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The luminosity is calculated using the following three 
equations relating apparent magnitude, flux and lumi- 
nosity at 4400 A (eg. Weedman 1998), 



F 



log ■ 



erg cm - ' s 



= -0.4(m s + 20.42) 

-i A 



L , =4itd 2 (l + z) 2 F s 

A, 4400 A v ' A, 4400 A 



c[go£ + (go - 1)(V1 + 2g a z - 1)] 
q 2 H (l + z) 



(3) 



(4) 



(5) 



Here, F s and L s are the median flux and lu- 

A, 4400 A A, 4400 A 

minosity at 4400 A, c is the speed of light in a vacuum. 
The apparent magnitudes, mj, are taken from G99 and 
corrected for the line-of-sight extinction, with B band ex- 
tinction values taken from Schlegel, Finkbeiner & Davis 
(1988). Here we use H = 70 km s" 1 Mpc" 1 and and 
go = 0.2 for consistency with the G99 preparation of the 
lightcurves. 

The corresponding errors in luminosity, disregarding 
the uncertainty in Ho, are typically ±1.6 x 10 41 erg s _1 
and are derived from the G99 estimate of a photometric 
accuracy for this band. 



2.1 Sample Refinement 

PG1226+023, also known as 3C 273, is the most optically 
luminous object in this sub-sample of the PG quasars, 
an order of magnitude greater in luminosity than all the 
other sources. This AGN is known to have optical, X-ray 
and radio jets. It is the only source in the optically se- 
lected sample which shows definite jets, though 7 of the 
42 objects are radio-loud. The jet emission is thought to 
be synchrotron and variable itself. Since we cannot sepa- 
rate the disc emission from that of the jet PG1226+023 
will not be included in this analysis. 
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3 THE MODEL 

We assume the optical emission observed in the PG sam- 
ple originates from an accretion disc. Adopting the ap- 
proach of SC89, the disc may be divided into two distinct 
regions: the optically thick outer region emitting black- 
body contributions into the B band, and the inner region • 
of the disc, which does not significantly contribute to the ' 
B band emission. 

We assume that the optical variability is caused by an in- 
stability on the thermal time-scale, which varies the po- 
sition of the inner radius of the optical emission region. 
We assume that the two-folding time-scales we measure 
in the B band lightcurves of the PG sample correspond 
to the thermal time-scale at the radius Rent where the 
luminosity in the B band has decreased by 50 per cent 
from its maximum value. The value of the radius -R cu t for 
a given luminosity depends upon the disc model chosen, 
accretion rate and black hole mass. 

Comparing -R cu t for each model with the thermal time- 
scale, Tth, gives a value of a from the following equation 

Tth = — (6) 
ato k 

where Wk = V^f? * s * ne Keplerian angular velocity. 

To describe the outer optical flux emitting disc we 
use a standard multi-temperature blackbody a-disc code 
developed by Czerny & Elvis (1987) and used in SC89. 
The model assumes a geometrically thin, optically thick 
Keplerian disc accreting steadily onto a Schwarzschild 
black hole. 

We assume that the disc radiates locally as a blackbody 
with an emitted flux of 
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where R is the radial distance from the central object, 
a is the Stefan-Boltzmann constant, T c ff is the effective 
temperature of the disc at radius 7?, G is the gravita- 
tional constant, M is the black hole mass, M is the 
steady accretion rate and -Rmner is the inner radius of 
the disc (eg. Frank, King & Raine 1992). We assume 
dinner = 3i? g , where R g = 2GM/c 2 . Inside the radius 
Rcut we do not make any further assumptions about the 
nature and spectral properties of the disc. We only as- 
sume that the optical flux from that region is negligible. 
Fig.^shows different spectra depending on the choice of 
Rcut- Note that the maximum contribution to the flux at 
a given frequency v comes from the disc region with the 
temperature kT c s — 1.65hv, which is located at the dis- 
tance R„ = 0.77xl0 24 //- 4/3 (M/M Q )- 1/3 (L/L E dd) 1/3 -R g 
(SC89), where L E dd ~ 1-3 x 10 38 (M/M Q ) ergs^ 1 is the 
Eddington luminosity. Hence for a source with a black 
hole mass of 10 s Mq emitting at the Eddington luminos- 
ity the location of the B band is R v ~ 21 — 38 R g , within 
our chosen range. 

The shape of the spectrum depends upon 

1) M, the black hole mass 

2) M, the accretion rate 

3) the extent of the outer disc contributing to the optical 
flux, which depends on the variable R cu t since -Router is 
fixed at 10 3 R a 




Figure 1. The model blackbody spectrum for a M = 10 8 Mq 
black hole, accreting at one tenth of the Eddington rate, with 
increasing radial disc size i? cu t = 3 (upper curve), 10, 29, 64, 
127, 218, 345, 514 (lower curve) R g . The outer disc radius is 
always fixed at -Router = 10 3 i? g . 

The vertical dotted lines mark the region covered by the B 
band in the observers' frame (3600-5600 A). 



4) i, the disc inclination. We assume an inclination angle 
of cosi = 0.75 (41°) since optical spectroscopy has shown 
that all the sample objects have broad emission lines. 

We also consider the effects of opacity in a modified 
blackbody model in which electron scattering in the disc 
is taken into account (Czerny & Elvis 1987; Haardt & 
Maraschi 1993). Electron scattering becomes important 
in the disc at high temperatures, so we expect the ef- 
fects of opacity to be visible at ultraviolet energies and 
above, but should not significantly modify the continuum 
emission at optical wavelengths. 

To enable comparison with the data the luminosities 
of the model spectra were convolved with the B filter 
which had been applied during observation of the PG 
sample. The integrated luminosity over the Johnson B 
band, Lb, is calculated and directly compared with the 
observed B luminosities. Two-folding time-scales for the 
models are readily obtained from the logL vs. logr t h 
plots by measuring the time taken for the luminosity 
to drop from its maximum to half that value (Fig. 
Each value of a for a set of models has a characteristic 
time-scale (Equation 6) . Finding sets of models which en- 
compass the data sample allows the determination of the 
range of the viscosity parameter for that data set. 



4 RESULTS 
4.1 Variability 

The lightcurves for this dataset (G99) show a typical in- 
trinsic variability amplitude of as — 0.14 mag with in- 
trinsic rms amplitudes of 5% < ob/B <34% over the 
seven year period. The variability time-scales lie approx- 
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Figure 2. This figure shows the blackbody disc models' B 
band luminosity vs. thermal time-scale. 

This is the case of L/L Edd = 0.1 and black hole masses 
10 7 , 10 75 and 10 s M , with the smallest mass producing the 
lowest luminosity. 

For each mass the model is calculated for 8 steps in R C ut from 
the last stable orbit (3-R g ) out to 512R g . This has been done 
for two values of the viscosity parameter; a = 0.01 (dashed 
lines) and a = 0.02 (dotted lines). The solid lines join the 
points at which the luminosity has decreased by 50 per cent 
from its maximum value in each model (t2). 



Figure 4. Comparison of the observed distribution of two- 
folding time-scales for 41 PG sample AGN (binsize A log T2 = 
0.2, N counts per bin) with the distributions predicted from 
random realisations of the broken power law power spectrum 
(with break frequency Vbk) which probably describes the red- 
noise optical variability in AGN (10000 lightcurves simulated, 
N rescaled to 41 counts over whole distribution). 
Solid histogram: observed distribution, dotted line: i/^k = 
10~ 7 Hz, dot-dashed line: !/bk = 2 X 10 -8 Hz, dashed line: 
u hk = 4 X 10" 9 Hz. 
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Figure 3. The two-folding time-scales for each source versus 
luminosity in the B band (log scale). The mean time-scale and 
standard deviation of the mean are represented by the solid 
line and dot-dashed lines respectively, with the vertical error 
bar at the point of the mean luminosity. The uncertainty in 
the value of Ho introduces uncertainties in the luminosities, 
so we have not calculated the error on the mean luminosity 
here. The length of the luminosity error bar includes 95 per 
cent of the data (2a). 



imately between 300 and 4500 days. The distribution of 
two- folding time-scales is shown in Figs. Inland 0] and the 
individual time-scales and luminosities for each AGN are 
given in Table 2. The mean two-folding time-scale for the 
sample is 1566.4 days and the corresponding standard 
deviation of the mean is 143.6 days (in log space, mean 
log t 2 = 3.12 ±0.04). 



4.2 Distribution of two-folding time-scales 

The optical variability of AGN is consistent with a red- 
noise process, so that each lightcurve is a realisation of 
an underlying process described by a power law power 
spectrum with a slope and amplitude determined by the 
mechanism which causes the variability. In principle, it 
is possible to derive the power-spectral parameters from 
the model for variability, and compare these with ob- 
servations, however for simplicity and ease of interpreta- 
tion we consider only the two-folding time-scale, which 
is a property of the power-spectral shape. Due to the 
stochastic nature of red-noise processes, lightcurves ob- 
tained at different times will look different, and show dif- 
ferent observed two-folding time-scales. Ideally, we need 
to measure the 'true' average two- folding time-scale by 
averaging the two-folding time-scale over many 7-year 
lightcurves from the same source. As this is clearly not 
possible, we must take account of the fact that the two- 
folding time-scale measured for a given source is not the 
same as the average two-folding time-scale, and consider- 
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Table 1. Individual minimum two-folding time-scales and lu- 
minosities for AGN in this PG subsample. The luminosity 
quoted is the median B band luminosity in the rest-frame 
of the source. Errors in the luminosities stem from errors in 
the observed magnitudes and hence the flux. Errors due to the 
uncertainty in Hq are not included here. 



Table 2. Model parameters for the PG data, with correspond- 
ing mass and accretion rate ranges of 10 7 < M < 10 8 ' 5 ( Mq) 
and 0.02 < M < 0.70 ( M year" 1 ). 



Quasar 


other names 


z 


RL 


L 4400A Xl ° 42 
(org s 1 ) 


T"2 

(days) 


PG0026+129 




0.142 




11.9±0.22 


1454 


PG0052+251 




0.155 




20.7±0.38 


318 


PG0804+762 




0.100 




14.7±0.27 


1009 


PG0838+770 




0.131 




1.33±0.02 


1907 


PG0844+349 


Ton 951 


0.064 




5.29±0.10 


567 


PG0923+201 


Ton 1057 


0.190 




11.1±0.20 


1158 


PG0953+415 


K 348-7 


0.239 




15.4±0.29 


1266 


PG1001+054 




0.161 




3.62±0.07 


301 


PG1012+008 




0.185 




16.6±0.31 


1994 


PG1048+342 




0.167 




1.29±0.02 


496 


PG1100+772 


3C 249.1 


0.313 


Y 


14.5±0.27 


3186 


PG1114+445 




0.144 




6.70±0.12 


1901 


PG1115+407 




0.154 




3.80±0.07 


2262 


PG1121+423 




0.234 




6.23±0.11 


793 


PG1151+117 




0.176 




6.64±0.12 


608 


PG1202+281 


GQ Com 


0.165 




5.02±0.09 


685 


PG1211+143 




0.085 




13.5±0.25 


1602 


PG1229+204 


Ton 1542 


0.064 




3.52±0.06 


1189 


PG1307+086 




0.155 




8.54±0.16 


904 


PG1309+355 


Ton 1565 


0.184 


Y 


10.4±0.19 


1874 


PG1322+659 




0.168 




6.10±0.11 


4329 


PG1351+640 




0.087 




9.92±0.18 


1891 


PG1354+213 




0.300 




3.69±0.07 


1953 


PG1402+261 


Ton 182 


0.164 




8.47±0.16 


2600 


PG1404+226 




0.098 




2.16±0.04 


1562 


PG1411+442 


PB 1732 


0.089 




7.35±0.14 


2750 


PG1415+451 




0.114 




2.90±0.54 


1116 


PG1426+015 


Mkn 1383 


0.086 




4.08±0.07 


1327 


PG1427+480 




0.221 




3.95±0.07 


1542 


PG1444+407 




0.267 




2.80±0.05 


3362 


PG1512+370 


4C 37.43 


0.371 


Y 


10.7±0.2O 


3194 


PG1519+226 




0.137 




3.38±0.06 


1566 


PG1545+210 


3C 323.1 


0.266 


Y 


10.3±0.19 


1470 


PG1613+658 


Mkn 876 


0.129 




12.9±0.24 


856 


PG1617+175 


Mkn 877 


0.114 




6.91±0.13 


356 


PG1626+554 




0.133 




5.91±0.11 


668 


PG1700+518 




0.292 




21.3±0.40 


985 


PG1704+608 


3C 351 


0.371 


Y 


22.5±0.41 


954 


PG2130+099 


II Zw 136 


0.061 




7.50±0.14 


1454 


PG2233+134 




0.325 




11.8±0.22 


2551 


PG2251+113 


PKS 2251+11 


0.323 


Y 


17.00±0.31 


2264 



able scatter is introduced into the distribution of observed 
two-folding time-scales for sources which have identi- 
cal variability properties (i.e. identical power-spectral 
shapes) . 

We can test whether the observed distribution of two- 
folding time-scales is consistent with a single power- 
spectral shape and amplitude (and a corresponding single 
average two- folding time-scale), by comparing the distri- 
bution we observe with that determined from simulated 
red-noise lightcurves generated from plausible power- 
spectral parameters. We therefore simulated random 
realisations of red-noise lightcurves using the method 
of Timmer & Konig (1995). To date, only the X-ray 
power-sp ectral shapes of AGN have been well con- 
strained ftjttlev. M c Hardv fc Papadakisl EoO^) . showing 
fairly steep (l/v 2 ) power-law slopes at high frequencies, 
breaking to flatter l/v slopes below some break fre- 
quency. The optical power-spectral shape of the AGN 
in our sample is not well constrained, but correlated 
long-time-scale optical/X-ray variations in the luminous 
Seyfert galaxy NGC 5548 suggest that although the opti- 



Model 


TIT 

^/^Edd 


<^min 


GJmean 


Omax 


disc blackbody 


0.01 


0.010 


0.011 


0.012 


disc blackbody 


0.1 


0.015 


0.017 


0.018 


disc blackbody 


0.25 


0.017 


0.019 


0.021 


disc blackbody 


0.5 


0.020 


0.022 


0.025 


disc blackbody 


0.75 


0.022 


0.024 


0.027 


disc blackbody 


1.0 


0.023 


0.026 


0.028 


modified blackbody 


0.01 


0.010 


0.011 


0.012 


modified blackbody 


0.1 


0.015 


0.017 


0.018 


modified blackbody 


0.25 


0.016 


0.018 


0.020 


modified blackbody 


0.5 


0.019 


0.021 


0.023 


modified blackbody 


0.75 


0.020 


0.022 


0.024 


modified blackbody 


1.0 


0.021 


0.023 


0.026 



cal power spectrum of luminous AGN is steeper than the 
X-ray power spectrum at hi gh frequencies . the sh apes are 
similar at low frequencies (lUttlev et al. 1 120031) . There- 
fore we assume high-frequency power-spectral slopes of 
l/v 2 ' 5 , breaking to l/v below some break frequency, f D k, 
which we assume is common to both X-ray and opti- 
cal power spectra. Breaks in the X-ray power spectra 
of AGN appear to scale linearly wi th black hole mass 
dUttlev. M c Hardv fc Papadakisll2002t . so we first chose a 
break frequency of 2 x 10 -8 Hz, corresponding to a black 
hole mass of a few x 10 8 Mq . We chose a power-spectral 
normalisation which corresponds to the typical observed 
~ 10 per cent variability amplitude in the sample. 
We simulated 10000 lightcurves of 7 year duration, and 
measured the resulting distribution of two-folding time- 
scale i~2, which is plotted for comparison with the ob- 
served distribution of T2 in Fig. [I] We also show the dis- 
tributions corresponding to break frequencies of 10~ 7 Hz 
and 4 x 10 -9 Hz, i.e. a possible factor ~ 25 range in 
black hole mass. The observed distribution of t~2 is con- 
sistent with that expected from random realisations of a 
single power spectral shape expected from few x 10 s M 
black holes. So the distribution of black hole masses is 
likely to be narrow (< decade) and the variability in the 
lightcurves of all the AGN can be characterised by a sin- 
gle identical two-folding time-scale, as expected from a 
homogeneous sample of AGN such as that presented here. 
This result is not sensitive to the slope above the break. 
From Equation 6 the sample must also be consistent with 
a single value of the quantity R^ 2 /a. Keeping the Ed- 
dington ratio L/L^dd fixed, there is a single value of the 
viscosity parameter a determined by the average T2 of 
the sample. 

For 0.01 < L/L-Edd < 1.0 the data are best fitted 
by the set of models outlined in Table 2. Fig. shows 
the allowed values of a over a range of the Eddington 
ratio for both the disc models we apply. We note that the 
difference in viscosity between blackbody and modified 
blackbody models is small over this wavelength range. 
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0.030 




0.030 




Figure 5. This figure shows the allowed range of a values for both the blackbody models (left) and the modified blackbody 
models (right). The permitted values of a lie within the solid lines, the best fitting values lying along the dot-dashed line. The 
shaded areas are regions in which a given viscosity parameter, a, and Eddington ratio, L/L Edd , cannot occur simultaneously in 
the model. 



5 DISCUSSION 

Since direct observations cannot yet resolve the inner ac- 
cretion disc in AGN, numerical simulations of accretion 
flows have existed for several decades and are moving in 
the direction of global 3-D models. Our results fall within 
the range 0.005 < a < 0.6 predicted by a variety of recent 
numerical simulations of MHD turbulent discs. However, 
we obtain much tighter constraints towards the lower end 
of this broad range which could eliminate some models. 
For example, Armitage (1998) writes that a ~ 0.1 in 
global and local cylindrical disc simulations with net ver- 
tical field, whereas a ~ 0.01 when the net vertical field 
is zero. MHD simulations have also shown that an al- 
ternative Q-prescription in which stress is proportional 
only to the magnetic pressure may be a better diagnos- 
tic of the viscosity mechanism, since using the Shakura 
Sunyaev (1973) a-prescription a rang e of a values can 
be ob tained for the same gas pressure dHawlev fc Balbusl 
1995). 



The only previous observational viscosity estimates 
for AGN discs were those calculated by SC89 for 8 AGN 
and 4 blazars at 1060 A , 3 AGN at 1395 A and 12 AGN 
at 1740 A , with the prediction that a m i n = 0.001 and 
awx ~ 0.1. Our results allow a much tighter estimate of 
the viscosity range and we find a value around 0.02. Both 
estimates have used the same disc model, but at different 
wavelengths; we chose an optical sample rather than the 
UV data used by SC89 since the optical bands lie further 
from the peak in the power law region of the disc spec- 
trum. SC89 assume an inclination angle of 0° and include 
blazars in their sample. They discuss the likelihood that 
the interpretation of blazar variability as thermal changes 
in the accretion disc may not be valid since synchrotron 
radiation is known to dominate in many of these objects. 
They also note that the sampling of 5 of their sample 



objects was not designed for detecting the time-scales of 
interest when estimating disc viscosities. O'Brien et al. 
(1988) examined IUE data of several AGN and found that 
most had UV two-folding time-scales of weeks to years 
with no shorter time-scale variability. The sampling rate 
for this data set was poor. Consequently, it would not be 
possible to estimate a for their sample. The PG subsam- 
ple used here is larger and has a higher and more even 
sampling than any sample used in previous studies of this 
kind. 

We find this sample of 41 AGN all have similar (< 
decade) masses and therefore can all be represented by a 
single two- folding time-scale. Hence we find a very nar- 
row range of a and/or L/L Ed d- This seems to suggest 
that we are observing all these AGN in a single 'state'. If 
this is true, it may be that all observed AGN are in the 
outburst state and those in the quiescent state are either 
too faint to have yet been detected or appear as nor- 
mal galaxies (Siemiginowska & Elvis 1997; Burderi, King 
& Szuszkiewicz 1998; Hatziminaoglou, Siemiginowska & 
Elvis 2001). If the former is true, this would then imply 
that the total number of AGN in existence is far higher 
than current estimates based on the number of sources 
we see. However, a single value of a throughout the ther- 
mal limit cycle has been predicted in some numerical 
simulations of MHD turb ulence in AGN accretion disks 
iMenou fc Quataertfe OOl). though the thermal ionisation 
instability is unlikely to be a source of large amplitude 
variability in AGN. 

Reverberatio n mapping measur ements exist for 16 of 
these sources dKaspi et al. 12 000) ranging from 2.16 x 10 
to 4.7 x 10 s M (M aV erage = 1.51 x 10 8 M Q ) with typi- 
cal uncertainties of about 40 per cent. Black hole masses 
we find through the disc modelling method are depen- 
dent upon the ratio of AGN luminosity to the Eddington 
Limit, L/LEdd. For L/L Edd = 0.01, M = 10 s - 10 8 5 M Q 
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as suggested by comparing the observed distribution of 
two-folding time-scales with distributions predicted from 
random realisations of the broken power law power spec- 
tra. The reverberation mapping results are compatible 
with all the Eddington ratios used here. 

We now summarise the main assumptions and note 
some caveats on our determination of a. The accretion 
disc models we use are simple multi-temperature black- 
body spectra (Czerny & Elvis 1987) , and we assume that 
the disc is the source of the B Band luminosity. We also 
assume that long term variability in this band occurs on 
the thermal time-scale at a radius corresponding to a 50 
per cent change in B Band luminosity. The limits on a 
are derived only for the sample as a whole since the vari- 
ability time-scale and mass of each individual object is 
not precisely known. The results are dependent on the 
validity of the assumptions stated above and require a 
single viscosity value throughout the outer disc. In prac- 
tice, viscosity is likely to be both radially and time depen- 
dent and variations in luminosity may occur on various 
time-scales and radii. For example, the a value we de- 
rive here would be a lower limit on the true value if the 
observed variability time-scales are not equal to but are 
longer than the thermal time-scale. 

Blackbody and modified blackbody models gave very 
similar (and, at accretion rates one tenth of Eddington 
or lower, the same) values of a. The small difference be- 
tween the blackbody models and the addition of opacity 
effects to these models shows that in the optical B band, 
away from the peak of the emitted flux, electron scatter- 
ing has only a small effect on the spectrum as expected. 
The disc models do not include warps, dumpiness or 
flares. Contributions to the variability from occurrences 
of flares is, however, likely to be on time-scales shorter 
even than the sampling interval of the lightcurves. 
G99 estimate pollution by the host galaxies to be ap- 
proximately 5 per ce nt based on observations with HST 
flBahcall et al. 1 11997), The emission line contribution to 
the optical flux is estimated in the same paper at 5-10 
per cent. The models applied here are for AGN contin- 
uum flux only and we have not reduced the observed flux 
to account for these since a 10 per cent luminosity in- 
crease will not greatly affect the viscosity estimates. We 
also note that our luminosity estimates are affected by 
the choice of Ho which, if incorrect, would introduce an 
error into our results. 

These models do not incorporate a corona, transition re- 
gion, self-gravitation or irradiation of the disc material. 
Coronae are thought to be the source of hard X-rays but 
could alter the shape of the optical spectrum as we may 
see reprocessing and/or scattering, and this may be im- 
portant in some AGN. The role of the corona in form- 
ing the continuum emission is discussed in Kurpiewski, 
Kuraszkiewicz & Czerny (1997) and a corona could be in- 
corporated in future determinations of the a-parameter. 

There are currently no conclusive observational or 
theoretical bounds on the a-parameter. We have carried 
out this study in the hope of providing another clue to 
the nature of viscosity in AGN. 



6 CONCLUSIONS 

We assume that optical variability of AGN observed on 
time-scales of months to years is caused by local insta- 
bilities occurring on the thermal time-scale in the in- 
ner accretion disc. The disc is modelled using multi- 
temperature blackbody models. In this case we constrain 
the viscosity to 

0.01 < a < 0.03 for 0.01 < L/L Edd < 1.0. 

We show this sample is consistent with a single a value 
at a given Eddington ratio, determined by the average 
characteristic time-scale of all 41 AGN. The range of a 
values we find lies within the range predicted by current 
numerical simulations of MHD turbulent discs. The mass 
range we obtain for this sample is consistent with the re- 
verberation mapping measurements made for 16 of these 
PG quasars. However, we do stress that this is a simpli- 
fied model of an accretion disc, and merely a first step in 
observational viscosity estimates for AGN. 
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